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Abstract 

We present a preliminary measurement of CP-violation parameters in the decay — > K~ , 
using approximately 347 million BB events collected by the BABAR detector at SLAC. Reconstruct- 
ing the neutral kaon as — > 7r"'"7r~,i^g — > vr'^Tr^, or K^, we analyze the Dalitz plot distribution 
and measure fractions to intermediate states. We extract CP parameters from the asymmetries 
in amplitudes and phases between B^ and B^ decays across the Dalitz plot. For decays to ipK^, 
we find f3eff = 0.06 ± 0.16 ± 0.05, Acp = -0.18 ± 0.20 ± 0.10, where the first uncertainty is sta- 
tistical and the second one is systematic. For decays to /oK°, we find Pejj = 0.18 ± 0.19 ± 0.04, 
Acp = 0.45 lb 0.28 ± 0.10. Combining all K~^K~K^ events and taking account of the different CP 
eigenvalues of the individual Dalitz plot components, we find (3eff = 0.361 it 0.079 it 0.037, Aqp = 
—0.034 lb 0.079 lb 0.025. The trigonometric reflection at it/2 — (S^jj is disfavored at 4.6(7. We also 
study angular distributions in B^ —>■ K~^K~Kg and B^ — > (pK^ decays and measure the direct 
CP asymmetry in B+ (j)K+ decays, Acp = 0.046 zb 0.046 zb 0.017. 
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1 INTRODUCTION 



Wc describe a i?-flavor tagged, time-dependent Dalitz plot analysis of the K^K~K^ de- 

cay [1], with the reconstructed as vr^vr^, vr'^vr'^, or K^. In the Standard Model 

(SM), these decays are dominated by 6 — > sss gluonic penguin amplitudes, with a single weak 
phase. Contributions from h uqq tree amplitudes, proportional to the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element Vub with a CP-violating weak phase 7 [2], are small, but may 
depend on the position in the Dalitz plot. In (j){K^ K~)K^ decays the modification of the 

CP asymmetry due to the presence of suppressed tree amplitudes is at 0(0.01) [3, 4], while at higher 
K'^K~ masses a larger contribution at O(O.l) is possible [5]. Therefore, to very good precision, we 
also expect the direct CP asymmetry for these decays to be small in the SM. The CP asymmetry 
in B'^ K^K~K^ decay arises from the interference of decays and B^ ^ B^ mixing, with a 
relative phase of 2/3. The Unitarity Triangle angle (3 has been measured in B^ \c(i\K^ decays to 
be sin 2/3 = 0.685 it 0.032 [6, 7]. Current direct measurements favor the solution of /3 = 0.37 over 
/? = 1.20 at the 98.3% C.L. [8]. 

The decay B^ K^K^K^ is one of the most promising processes with which to search for 
physics beyond the SM. Since the leading amplitudes enter only at the one-loop level, additional 
contributions from heavy non-SM particles may be of comparable size. If the amplitude from heavy 
particles has a CP- violating phase, the measured CP- violation parameters may differ from those 
expected in the SM. 

Previous measurements of the CP asymmetry in B^ K^K^K^ decays have been performed 
separately around the (f) mass, and for higher K'^K~ masses, neglecting interference effects between 
intermediate states [9]. In this analysis, we extract the CP-violation parameters by taking into 
account different amplitudes and phases across the B^ and B^ Dalitz plots. 

2 EVENT RECONSTRUCTION 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy B factory at SLAC. A total of 347 million BB pairs were used. 

The BABAR detector is described in detail elsewhere [10]. Charged particle (track) momenta 
are measured with a 5-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) coaxial with a 1.5-T superconducting solcnoidal magnet. Neutral cluster (photon) positions 
and energies are measured with an electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) 
crystals. Charged hadrons are identified with a detector of internally reflected Cherenkov light 
(DIRC) and specific ionization measurements (dE/dx) in the tracking detectors (DCH, SVT). 
Neutral hadrons that do not interact in the EMC are identified with detectors, up to 15 layers 
deep, in the flux return steel (IFR). 

We reconstruct B^ —>■ K~^K~ decays by combining two oppositely charged tracks with a 
ifg — TT+TT", — 7r°7r°, or candidate. — ^ (f)K'^ decays are reconstructed from three 
charged tracks. The and K'^ tracks must have at least 12 measured DCH coordinates, a 
minimum transverse momentum of 0.1 GeV/c, and must originate from the nominal beam spot. 
Tracks are identified as kaons using a likelihood ratio that combines dE/dx measured in the SVT 
and DCH with the Cherenkov angle and number of photons measured in the DIRC. For K'^K~ 
masses near the mass, higher efficiency kaon identification is used, while for higher masses higher 
purity criteria are chosen to reduce background. 

For all modes, the main background is from random combinations of particles produced in 
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events of the type e+e" qq {q = u,d,s,c) (continuum). Additional background from decays of 
B mesons to other final states, with and without charm particles, is estimated from Monte Carlo 
simulation. 

We use event-shape variables, computed in the center-of-mass (CM) frame, to separate con- 
tinuum events with a jet-like topology from the more isotropic B decays. Continuum events are 
suppressed with a requirement on the quantity cos(^t), cos(^t) < 0.9, where is the angle be- 
tween the thrust axis of the B candidate's daughters and the thrust axis formed from the other 
charged and neutral particles in the event. We select events in the range . Further discrimination 
comes from the Legendre moments >Ci=o,2 = J2j PjLi{Oj), where the sum is over all tracks and 
clusters not used to reconstruct the B meson; Lj is the Legendre polynomial of order i, and 6j 
is the angle to the B thrust axis. Lastly, the magnitude of the cosine of the angle of the B with 
respect to the collision axis | cos6b\, is also used. 

2.1 B° K+K K^, TT+TT- 

For decays B^ K^K^Kg and Kg — > tt+tt^, candidates are formed from oppositely charged 
tracks with an invariant mass within 20 McV/c^ of the K° mass [2]. The Kg vertex is required to 
be separated from the B^ vertex by at least 3a. The angle a between the Kg momentum vector 
and the vector connecting the B^ and K^ vertices must satisfy cos a > 0.999. 

B candidates are identified using two kinematic variables that separate signal from continuum 
background. These are the beam-energy-substituted mass mss = \J {s/2 + Pj • pb)^/-E? — p^, 
where ^/s is the total e+e" CM (CM) energy, (£'i,pi) is the four-momentum of the initial e+e" 
system and pe is the B candidate momentum, both measured in the laboratory frame, and AE' = 
Eb ^ \/s/2, where Eb is the B candidate energy in the CM frame. Distributions of these variables 
in data, for signal and background events calculated using the sVlot event- weighting technique [11], 
are shown in Fig. 1. 

2.2 B« K+K-K^, ttOtt" 

For decays B^ — K'^K~Kg and Kg — 7r°7r°, Kg candidates are formed from two 7r° — 77 
candidates. Each of the four photons must have > 0.05 GcV and have a transverse shower shape 
loosely consistent with an electromagnetic shower. Additionally, we require each vr^ candidate 
to satisfy 0.100 < m^^ < 0.155 GeV/c^. The resulting vr^TT^ mass is required to satisfy 

0.4776 < m^o^o < 0.5276 GeV/c^. A Kg mass constraint is then applied for the reconstruction of 
the B^ candidate. 

The kinematic variables mss and AE are formed for each candidate as in Sec. 2.1. Distributions 
of these variables are shown in Fig. 2. Note that the mean of the signal AE distribution is shifted 
from zero due to energy leakage in the EMC. 

2.3 B° K+K-K^ 

We identify a K^ candidate either as a cluster of energy deposited in the EMC or as a cluster 
of hits in two or more layers of the IFR that cannot be associated with any charged track in the 
event. The K^ energy is not measured. Therefore, we determine the K^ laboratory momentum 
from its flight direction as measured from the EMC or IFR cluster, and the constraint that the 
invariant K~^K~K'^ mass equal the B^ mass [2]. In those cases where the K^ is detected in both the 
IFR and EMC we use the angular information from the EMC, because it has higher precision. In 
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Figure 1: Distributions of kinematic variables (left) m-ES and (right) AE for the K~ Kg{7:^'K^) 
subsample: (top) signal, (bottom) continuum background. The points are data events weighted 
with the s'Plot technique [11], and the curves are the PDF shapes used in the ML fit (Sec. 3). 
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Figure 2: Distributions of kinematic variables (left) m^s and (right) AE for the K'^^tt^tt^) 
subsample: (top) signal, (bottom) continuum background. The points are data events weighted 
with the s'Plot technique [11], and the curves are the PDF shapes used in the ML fit (Sec. 3). 
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order to reduce background from ir^ decays, we reject an EMC candidate cluster if it forms an 
invariant mass between 100 and 150 MeV/c^ with any other neutral cluster in the event under the 
77 hypothesis, or if it has energy greater than 1 GeV and contains two shower maxima consistent 
with two photons from a vr" decay. The remaining background of K'^ candidates due to photons and 
overlapping showers is further reduced with the use of a selector based on the Boosted Decision Trees 
algorithm [12, 13]. This selector is constructed from cluster shape variables, trained with Monte 
Carlo events, and tested on Initial State Radiation e+e~ — (f){-^ KgK^)j events and reconstructed 
J/ipK^ candidates, which give a very pure sample of candidates. 

We use the kinematic variable A£' to characterize the signal and background. This variable, 
computed after the mass constraint on the B^ candidate, has a resolution of about 3.0 MeV for 
EMC events, and about 4.0 MeV for IFR events. For signal events, is expected to peak at 
zero, with a broad tail for positive values of /S.E. We require /S.E < 30 MeV, in order to be able 
to determine the shape of background under the signal peak. The mean value and the resolution 
of this variable has been taken from reconstructed B° — J/ipK'l events. 

In addition to the shape variables described in Sec. 2, we consider other variables, related to 
the missing energy in the event, to characterize events to the final state. The first, p^^^^^, is the 
difference between the T(45) energy and the total measured energy of the event, not including the 

candidate, projected onto the plane transverse to the beam axis. We also use the angle between 
p^jgg and the reconstructed direction, and the difference between the total visible energy of 
the event and the two reconstructed kaon energies. The latter corresponds to the unmeasured -fC^ 
energy, and has a somewhat diff'erent distribution in signal than in continuum background events. 
These three variables are used as inputs to a Fisher discriminant which has been trained on Monte 
Carlo samples and validated on data control samples. 

We optimized the selections on all of these variables for maximum sensitivity in the measurement 
of CP parameters. The selection is optimized independently in the region mj^+j^- < 1.1 GeV/c^ 
and in the rest of the Dalitz plot, because the higher K^K^ mass region gets more background 
from low momentum neutral candidates, for which the separation between photons and is worse. 
The final average eflficiency of the selection is about 25% in the low-mass region and 10% for the 
rest of events. 

The /S.E distribution, with the result of the fit superimposed, is shown in Fig. 3, after making 
a requirement on the ratio of signal likelihood to signal-plus-background likelihood to enhance the 
signal. 

2.4 B+ cf)K+ 

We use the B~^ (pK~^ decay to measure the charge asymmetry, the P-wave fraction, and the rela- 
tive phase between the S and P waves in the 0(1020) region. The selection of the (p meson candidate 
is done by applying an invariant K'^K~ mass cut defined as 1.0045 < rrix+K- < 1-0345 GeV/c^. 
For the bachelor candidate from the B~^ decay the track requirements are the same as for the 
(f> daughters but we apply a more restrictive kaon identification criterion. 

The distributions of mss and nix+K- variables after a cut on the ratio of signal and background 
probabilities is shown in Figure 4. 
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Figure 3: Distribution of the kinematic variable AE for the subsample. The sohd hne 

represents the total likelihood, while the dashed line represents the sum of continuum and BB 
background. A requirement on the ratio of signal likelihood to signal-plus-background likelihood 
is applied to enhance the signal, with an efficiency of about 30% for signal. 
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Figure 4: Distributions of the event variables (a) rriES and (b) mx+K- the (j)K~^ final state after 
reconstruction and a requirement on the likelihood calculated without the plotted variable. The 
efficiency for the selection and likelihood requirements is 78% for (a) and 95% for (b). The solid 
line represents the fit result for the total event yield and the dotted line for the background. 
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3 ANALYSIS OF THE DALITZ PLOT 



We analyze selected events in all neutral B samples using a maximum likelihood fit with the 
likelihood function £ for each subsample, defined as 



where j runs over all events in the sample, and rii is the yield for event category i. The probability 
density function (PDF) V is formed from multiple observables as 



Here cos Oh is the cosine of the helicity angle between the and the in the K^K~ center- 
of-mass frame, qtag is the flavor of the initial state, and At = tree — Uag is the difference of the 
proper decay times of the two B-mesons. a At is the error on At. The PDF for the time-dependent 
Dalitz plot, VdPi is described in detail later in the text. 7^ is a standard At resolution function 
with parameters evaluated in exclusive decays into final states with a charm meson as in our 
CP-asymmetry measurements in J/ipKg decays [6]. For the K'^K~K^ submode, the mES variable 
does not enter the likelihood function defined in Eq. (2). 

In all K~^K~K^ submodcs, the signal components of the PDFs for V(mEs) and V(AE) are 
parameterized using modified Gaussian distributions: V = exp[— (x — xo)'^/{2a'^ + a±{x — xq)^)], 
where x is the dependent variable. We determine the parameters xq, cr+, a-, a+, and a_ using 
simulated events, and fix them in fits to data. For x < xq {x > xq), the parameters (T_, a_ (cr+, a+) 
are used. In the B^ <pK~^ mode, the signal rriES distribution is parameterized as above, but 
with a± set to zero. The signal AE distribution is parameterized as the sum of two Gaussian 
distributions. 

In all applicable submodes, we use the Argus function to model the continuum background 
component of P(mEs) [14]. For P(A£'), linear shapes are used for the continuum background, 
except for the K^K^K'l submode. In that case, we use a reflection of the Argus function. 

In the Kg^TT^TT^) and K~ K'^ submodes, BB background components are parameter- 

ized with the same functional forms as the continuum backgrounds. Due to non-negligible corre- 
lation between twes and AE for BB background in the K~ Kg{'K^Tr'^) submode, we construct 
that PDF component as a two-dimensional histogram PDF in those variables. 

3.1 Background Decays in the Time-Dependent Dalitz Plot 

The Dalitz plot for the continuum background is parameterized using a two-dimensional histogram 
PDF in the variables mx+K- and cos 9h- The histogram is filled with candidates from the region 

5.2 < mES < 5.26 GeV/c^ for the K^K~K^ submodes. For the K^K~K^ submode, candidates 
from the region 20 < AE < 40 MeV are used. The At distribution is described with a separate PDF 
which, similarly to our previous measurement [9], uses a double-Gaussian resolution function and 
allows a fraction of decays to have a non-zero lifetime. For the S+ (f)K^ mode, the continuum 
backgroTind distribution in rrij^+j^- is modeled with the sum of a relativistic Breit-Wigner function 
(see Section 3.2) and a second-order polynomial. 

We estimate the amount of BB background from Monte Carlo events and again describe the 
Dalitz plot using a two-dimensional histogram PDF. The BB background in K'^K~Kg modes 
is almost purely combinatorial and is a few percent of the total background. The K'^K~K^ 




(1) 



V = r{mEs) • V{AE) • VDp{mK+K- , cos 9h, At, qtag) O ^(At, aAt)- 



(2) 
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sample, in addition to combinatorial BB background, contains decays into 4-body final states 
where a pion is missed in the reconstruction (Sec. 2.3). Such decays (e.g. K'^K~K*) are simulated 
with the assumption that the phase space contains the same distribution in the Dalitz plot as in 
B+ K+K+K- and B° K+K'K'^ decays: (pK*, a wide scalar at 1.5 GcV/c^ and a non-flat 
non-resonant distribution. The At distribution is described as a separate PDF that has a non-zero 
lifetime. The time-dependent CP asymmetry of this PDF, set to zero in the reference fit, is varied 
as a systematic uncertainty. 

The decays B^ D+R- {D+ K+K^) and 5° ^ R- {D+ K+K^) are kinematically 
indistinguishable from signal decays. Wc include non- interfering amplitudes for these modes in our 
Dalitz plot model, parameterizing the mesons on the Dalitz plot as Gaussian distributions 
with widths taken from studies of simulated events. 



3.2 Signal Decays in the Time-Dependent Dalitz Plot 

When the flavor of the initial state qtag^ and the difference of the proper decay times At, are 
measured, the time- and flavor-dependent decay rate over the Dalitz plot can be written as 



1 1 e~\^^\/'^BO 



1^1^ + 1^1^ ± qtag 2Im {AA*) sin AmrfAt (3) 

- qtag (l^P - 1^1^) COS Am^At ] , 

where plus (minus) sign is for decays to K+R-R^ {R+R^RD and qtag = +1 (-1) when the other 
B meson is identified as a SO (SO) using a neural network technique [6] . Approximately 75% of the 
signal events have tagging information and contribute to the measurement of CP violation param- 
eters. After accounting for the mistag rate, the effective tagging efficiency is (30.4 ±0.3)%. Events 
without tagging information are still included in the fit as they contribute to the determination of 
the Dalitz plot parameters. Decay amplitudes A and A are defined in (5) and (6) below. M^o, 
Tgo, and Am^ are the mass, lifetime, and mixing frequency of the B^ meson, respectively [2]. 

Four-momentum conservation in a three-body decay gives the relation M^o + ml + + 
m\ = + JTzfg -I- mla, where mfj = {pi + Pj)'^ is the square of the invariant mass of a 

daughter pair. This constraint leaves a choice of two independent Dalitz plot variables to describe 
the decay dynamics of a spin-zero particle. In this analysis wc choose the R^R^ invariant mass 
nix+K- cind the cosine of the helicity angle cos 9h- The PDF for the Dalitz plot rate becomes 

Vdp oc dViniK+K- , cos 0h, At, qtag) ■ £{mK+K- , cos Oh) • | J| (8) 7^(A^, a At), (4) 

where the Jacobian \ Jimj^+j^-)\ = {2mj^+j^-)(2qp) is given in terms of the charged kaon momen- 
tum q and neutral kaon momentum p, in the R^R~ frame. The efficiency e is calculated from 
high-statistics samples of simulated events and depends on the position on the Dalitz plot. 

The amplitude A {A) for the decay S° R+R'R^ (S° R'R+R^) is, in our isobar model, 
written as a sum of decays through intermediate resonances: 

A = 5^cv(l+6^)e^(<^'-+'''-+^) -Z^, and (5) 

r 

A = ^c,(l-6,)e^('^^-^'-'^)-/,. (6) 

r 

The parameters and are the magnitude and phase of the amplitude of component r, and 
we allow for different isobar coefficients for B^ and B^ decays through the asymmetry parameters 
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br and 6r- The parameter /3 is the CKM angle /3, coming from B^-B^ mixing. The function 
fj. = FrXTrX Zr describes the dynamic properties of a resonance r, where Fr is the form-factor 
for the resonance decay vertex, is the resonant mass-Uneshape, and Zr describes the angular 
distribution in the decay [15, 16]. 

Our model includes the ^(1020), where we use the Blatt-Weisskopf centrifugal barrier factor 
Fr = l/-\/l + {Rq)^ [15], where q is the daughter momentum in the resonance frame, and R is the 
effective meson radius, taken to be i? = 1.5 GeV (0.3 fm). For the scalar decays included in our 
model (/o(980), Xo(1550), and Xco)i wc use a constant form-factor. Note that we have omitted a 
similar centrifugal factor for the B^ decay vertex into the (pK^ intermediate state since its effect is 
negligible due to the small width of the c/){1020) resonance. 

The angular distribution is constant for scalar decays, whereas for vector decays Z = —Aq ■ p, 
where q is the momentum of the resonant daughter, and p is the momentum of the third particle 
in the resonance frame. We describe the line-shape for the 0(1020), Xo(1550), and XcO using the 
relativistic Breit-Wigner function 

T{m) = — 2 ^ — FT^' (7) 

nir — m|^+^_ — iniri. [m) 

where rrir is the resonance pole mass. The mass-dependent width is given as T{mx+K-) — 
Tr {q/qr)'^^~^^ (iT^r/^K+K-) (-^r (5)/-^r (<?r))^ ) where L is the resonance spin and q = qr when nij^+x- 
rrir. For the 0(1020) and XcO parameters, we use average measurements [2]. The Xo(1550) is less 
well-established. Previous Dalitz plot analyses of B+ K+K+K' [17, 18] and 5° K+R-R^ 
decays [19] report observations of a scalar resonance at around 1.5 GeV/c^. The scalar nature has 
been confirmed by partial- wave analyses [9, 18]. However, previous measurements report inconsis- 
tent resonant widths: 0.145±0.029 GeV/c^ [17] and 0.257 ±0.033 GeV/c^ [18]. Branching fractions 
also disagree, so the nature of this component is still unclear [20]. In our reference fit, we take the 
resonance parameters from Ref. [18], which is based on a larger sample of BE decays than Ref. [17], 
and consider the narrower width given in the latter in the systematic error studies. 
The /o(980) resonance is described with the coupled-channel (Flatte) function 

T{mK+K-) = 2 \ V' (8) 



where PK{mK+K-) = 2^1 - 4m|-/m|-+^_, p^{mK+K') = 1 - 4m2/m|-+^_, and the coupling 
strengths for the RR and tttt channels are taken as g.,^ = 0.165±0.018 GeV/c^, qk/q-k = 4.21zb0.33, 
and rur = 0.965 ± 0.010 GeV/c^ [21]. 

In addition to resonant decays, we include non-resonant amplitudes. Existing models consider 
contributions from contact terms or higher-resonance tails [22, 23, 5], but they do not capture 
features observed in data. We rely on a phenomenological parameterization [17] and describe the 
non-resonant terms as 

Anr = (ci2e^'^i2e-°^'"?2 + ciae'-^^e-amfg ^ C2^e''t'^^ e''^^^) ■ (1 + 6jvii) • e^(^+'^~«) (9) 

and similarly for Anr- The slope of the exponential function is consistent among previous mea- 
surements in both neutral and charged B decays into three kaons [17, 18, 19], and we use a = 
0.14 ±0.02 GeV-2-c^. 



16 



We compute the direct CP-asymmetry parameters for resonance r from the asymmetries in 
amphtudes (6^) and phases {6r) given in Eqs. (5, 6). We define the rate asymmetry as 



and Peff{f) = P + dr is defined as the total phase asymmetry. The fraction for resonance r is 
computed 

^ Jdcos9H dniK+K- ■ \J\ ■ (lA-P + .... 
Jdcos9H dmK+K--\J\-{\A\^ + \A\^) ' ^ ^ 

The sum of the fractions can differ from unity due to interference between the isobars. 

3.3 Calculation of Angular Moments 

As an alternative description that is less dependent on a resonance model, wc analyze the Dalitz 
plot in terms of moments of the cosine of the helicity angle, cos 6h, in -B"*" — > (pK~^ and — > 
K'^K~Kg{TT~^TT~) decays. We only assume that the total amplitude is the sum of the two lowest 
partial waves, and ignore direct CP- violation effects: 

A (A) = AsPo{cos0h) ± e''t'ApPi{coseH), (12) 

where (j) = (j)p — (f)s is the relative phase between S and P-wave strengths. The it sign corresponds 
to or B^ decays, respectively. Pq and Pi are Legendre polynomials that describe the amplitude 
dependence on cos Oh for S'-wave and P-wave decays, respectively. Integrating Eq. (4) over At, the 
tag-dependent decay rate can be written in terms of Legendre polynomials as follows: 



dr 



dcosdn ■ drriK+K- " |«^| 



^ {Pi) X Piicosen) (13) 

=0,1,2 

xPo{coseH) + ^J^Ajx P2{cos en) 

cos(t> X Pi{cos6h), (14) 



qtag ■ {D) ^ "^AsA^ 
(AmdTBo)2 + 1 ^ V2 

where the flavor tagging is necessary to measure the S-P wave interference term proportional to Pi . 
In the decay B^ , this term is diluted by the imperfect tagging {D) and the B^ — B^ 

mixing. In charged B decays to the K~^K^K^ final state we set Arrid = and {D) = 1 since 
qtag corresponds to the charge of the final state. Using our data sample, we compute the Legendre 
moments {Pi), 

{Pi) ^ ^PKcos0^^,i) >V(i)/e(i), (15) 

i 

where W is the weight for event i to belong to the signal sample [11]. These weights are computed 
from maximum likelihood fits that do not use the mass or the helicity angle in the fit. Finally, the 
fraction of P-wave decays is computed as fp = v'5/4 • (P2) / (Pq) [9]. 
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4 RESULTS 



4.1 Dalitz Plot and Angular Moments 

In order to determine parameters of the Dalitz plot model, we perform a fit to 3091 

K~ K'^{7r^'n'^) candidates in the full Dalitz plot. In this step we assume that all decays have 
the same CP-asymmetry parameters. We vary the event yields, isobar coefficients of the Dalitz plot 
model, and two CP-asymmetry parameters averaged over the Dalitz plot. We find a signal yield 
of 879 lb 36 events. The isobar amplitudes, phases, and fractions are listed in Table 1. The sum of 
resonant fractions in our DP model is different from 100% due to interference between resonances. 

We compare our fractions with other Dalitz plot analyses using fiavor symmetry [24]. We 
find consistent fractions for decays through the (/'(1020) resonances with the — K'^K'^K~ 
decay [17, 18]. The fraction of /o(980)K° decays is consistent with our S+ — K'^K'^K" analysis, 
and all K~^tt~^it~ Dalitz plot analyses [17, 18, 25]. The fraction of non-rcsonant decays, 

which is predicted to be half of the contribution in B^ K~^K^K~ [24], is harder to compare 
since existing measurements in the charged mode are inconsistent. Our result agrees well with 
BABAWs result [18], and is within two standard deviations of Belle's result [17]. Determination of 
the wide scalar resonance at 1.5 GeV/cP, labeled as Xo(1550), is even more uncertain. Using the 
same resonant parameters we find a much smaller fraction than in BABAR's analysis [18], but our 
solution is more consistent with Belle's B^ — > K^K^K~ analysis [17]. 

In order to achieve better statistical precision in measurements of CP-asymmetry parameters 
that are described in the following sections, we combine the K~ Kg{Tr~^Tr~) sample with samples 
of 1599 B° ^ (7r°7r°) and 22341 B^ K+K-R^ candidates. Fixing the coefficients of 

the isobar model to the values extracted from the K^{'k^'k^) submode, we find signal yields 

of 138 ± 17 and 499 ± 52 events in the if+K-K°(7r°7r°) and K+K-K^ submodes, respectively. 
Projection plots for the Dalitz plot variables are shown for all three submodes in Figure 5. 



Decay 


Amplitude cv 


Phase 


Fraction JF^ (%) 


0(lO2O)ii:" 


0.0098 ± 0.0016 


-0.11 ±0.31 


12.9 ± 1.3 


/o(980)/CO 


0.528 ± 0.063 


-0.33 ±0.26 


22.3 ±8.9 


Xo(1550)i^o 


0.130 ± 0.025 


-0.54 ±0.24 


4.1 ± 1.8 


NR {K+K-) 


1 (fixed) 


(fixed) 






0.38 ±0.11 


2.01 ± 0.28 


91 ±19 


{K-K^) 


0.38 ±0.16 


-1.19 ±0.37 




XcoK^ 


0.0343 ± 0.0067 


1.29 ±0.41 


2.84 ±0.77 


D+K- 


1.18 ±0.24 




3.18 ±0.89 


D+K- 


0.85 ± 0.20 




1.72 ± 0.65 



Table 1: Isobar amplitudes, phases, and fractions from the fit to the B^ K~ Kg{7r~^7r~) 

sample. Three rows for non-resonant contribution correspond to coefficients of exponential functions 
in Eq. (9), while the fraction is given for the combined amplitude. Errors are statistical only. 

As an additional crosscheck of our Dalitz plot model, we compute angular moments and ex- 
tract strengths of the partial waves in K'^K~ mass bins using the B^ — <pK'^ and B^ — 
K'^ K~K g{Tr'^Tr~) samples. In this approach we rely only on the assumption that the two low- 
est partial waves are present, but make no other assumption on the decay model. Wc confirmed 
the non-existence of higher partial waves by determining that higher angular moments ({P3-5)) are 
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Figure 5: Distributions of the Dalitz plot variables (left) mj^+j^- and (right) cos 9h for signal 
events (points) compared with the fit PDF in the following sub-samples: a) K'^ Kg{iT^TT^), b) 
K~^K~ Kg{7r^7r^), c) K~ K^. All sub-samples use the same Dalitz plot model but have different 
efficiencies so the resulting distributions differ. 
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Figure 6: (pK^: (left) The relative P-wave fraction fp in the interval 1.0045 < mx+K- < 
1.0345 GeV/c^. (right) The moment (Pi) calculated with the ^(1020) helicity angle defined with 
respect to the kaon of the same charge as the B meson. The dashed line corresponds to the fit 
result. 



consistent with zero. 

Prom the fit to 4947 S+ — candidates, we find 624 ± 30 signal candidates in the mass 
region 1.0045 < rriK+K- < 1-0345 GeV/c^. The event weights Wj are computed from the likelihood 
without mj^+j^- and cos 9h- From (Pq) and (P2) we obtain the average fraction fp = 0.891 it 0.014. 
The distribution of fp in four bins of uix+k- is shown in Fig. 6. In order to determine the 
relative phase between S- and P- waves, we make a minimization of the moment (Pi) given with 
Eq. (14), by varying the phase while keeping S- and P-wave strengths fixed. For the /S-wave we 
assume negligible energy dependence in the (/)(1020) region and parameterize the resonance shape 
with a relativistic Breit-Wigner function. We estimate a phase difference between P- and S-wave 
of (78 lb 20)°; the systematic error due to the choice of the S'-wave model and the scan method 
is negligible. The moment (Pi) from data, with the best phase solution superimposed, is displayed 
in the right-hand plot of Figure 6. 

Similarly, we compute the P-wave strength as a function of K~^K~ invariant mass in 

po ^ 

K'^ K~ K^{iT^ 'K~) decays as shown in Figure 7. We find the total fraction of P-wave in our sample 
is 0.29 ± 0.03 (stat) integrated over the entire Dalitz plot, which is consistent with our previous 
measurement [9]. In our model, this P-wave contribution originates from (j){lQ2Q)K^ decays, and 
non-resonant events with K^K^ and K~Kg mass dependence that reflects into an effective P-wave 
that enhances the central part of the plot in Figure 7. 

4.2 CP Asymmetry in the how-K'^K~ Mass Region 

In order to measure CP-asymmctry parameters for components with low-K'^ K~ mass with reduced 
model-dependence from the rest of the Dalitz plot, we select events using a cut of ruj^+x- < 
1.1 GeV/c^. After this selection, there are 836 K+ K' K^iTr+Tr') candidates, 202 K+ K' K^in^Tr^) 
candidates, and 4923 K'^K~Kl candidates remaining. The most significant contribution in this 
region comes from (f)(1020)K^ and /o(980)i^'' decays, with a smaller contribution from a low- 
K~^K~ tail of non-resonant decays. We vary the isobar parameters for the 0(1020) and fix all other 
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Figure 7: K~ K^{'k^ 'k~): The absolute strength of P-wave decays as a function of K~^K~ mass. 
The points are signal-weighted data and the histogram corresponds to the Dalitz plot model. 



components to the results of the full Dalitz plot fit. We also vary the CP amplitude and phase 
asymmetries for the (/>(1020) and /o(980). The asymmetry for the other components is fixed to the 
SM expectation. We perform a fit to the K~ K^{t:^ ■k~) submode, then perform an additional 
fit to the entire K^K~K^ sample. We find signal yields of 252 it 19, 35 it 9, and 195 it 33 events for 
K^{t:^ t:^) , Kg{ir^7r'^), and K^K^K'^ respectively. Fig. 8 shows projections of the 

Dalitz plot distributions of events in this region. The CP-asymmetry results are listed in Table 2; 
the systematic uncertainties will be described in Sec. 5. The left plots in Fig. 9 show distributions of 
At for i3°-tagged and S^-tagged events, and the asymmetry A{At) = {N^o — N^o)/{Nqo + N^o), 
obtained with the s'Plot event- weighting technique [11]. Correlation coefficients r between CP 
parameters, found in the fit to the combined sample, are also shown in Table 2. 



Name 


Fitted Value 

K+K-K'^ 

(vr^vr^) Combined 


Correlation 
(Combined) 
12 3 4 


1 Acp{(t>K'') 

2 Peffm'') 

3 Acp{hK'') 

4 (3,ffifoK^) 


-0.10 ±0.23 -0.18 ±0.20 ±0.10 
0.02 ±0.16 0.06 ±0.16 ±0.05 
0.36 ±0.33 0.45 ±0.28 ±0.10 
0.04 ±0.18 0.18 ±0.19 ±0.04 


1.0 -0.08 -0.27 0.11 
1.0 0.46 0.74 
1.0 0.20 
1.0 



Table 2: CP-violation parameters for —>■ K^K^K^ for uix+k- < 1-1 GeV/c^. For the combined 
K'^K^K^ sample, the first error is statistical and the second is systematic. For the K'^K^Kg 
submode, only statistical errors are shown. 

The decay <t>K^i with highly suppressed tree amplitudes, is, in terms of theoretical 

uncertainty, the cleanest channel to interpret possible deviations of the CP-violation parameters 
from the SM expectations. Values of [5eff are consistent with the value found in [cc] decays [6, 7]. 

As a consistency check we compare Peff{<j)K^) against a quasi-two-body approach (Q2B) that 
selects events around the (h resonance, 1.0045 < rriK+K- < 1.0345 GeV/c^, in order to increase 
the purity of P-wave decays. We add separate cos 9h and rrix+K- PDFs to the likelihood to 
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Figure 8: For the \ow-K^K mass Dalitz plot fit, distributions of the DaUtz plot variables (left) 
rrix+x- and (right) cos 6h for signal events (points) compared with the fit PDF in the following 
sub-samples: a) K+ R- K^{tt+tt-), b) K+ R- K^{tt°tt'^), c) K+K-K^. 
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At (ps) At (ps) 

Figure 9: (top) At distributions and (bottom) asymmetries in the K~ Kg{TT^TT^) submode for 
(left) 1.0045 < mx+K~ < 1-0345 GeV/c^ and (right) the whole Dalitz plot. For the At distributions, 
B^- {B^-) tagged signal-weighted events are shown as filled (open) circles, with the PDF projection 
in solid blue (dashed red). 



further suppress the S'-wave decays. From this fit we find 166 it 15 (j)K^ and 151 it 22 <j)K^ signal 
events. The CP parameters in the channel B^ — > cpK^ only are: = 0.10 it 0.29 and C^k = 

0.28 ± 0.20; in the channel S° ^ S^k = 0.69 ± 0.35 and C^x = -0.28 ± 0.33, with 

statistical errors only. The simultaneous quasi-two-body fit to both (pK^ and flavor decay modes 
yields the result S'^i^ = 0.39 ± 0.23 and CrpK = 0.10 ± 0.18 which is fully consistent with the 
Dalitz plot result. In this comparison we neglect interference effects and use the approximation 
ScfiK ~ sin(2/3ej(f )(1 — &^)/(l+&^) and C^k ~ —Acp to relate the Dalitz plot CP-violation parameters 
to the Q2B CP-violation parameters for the <pK^ decay. 

We also measure the direct CP asymmetry in B^ (j)K~^ decays, defined as Acp = {^b- ~ 
Nb+)/{Nb- + In a fit for the <pK+ and (pK' yields we find Acp = 0.046 ± 0.046 ± 0.017. 

4.3 Average CP Asymmetry in B° K^K~K^ 

We fit the average CP-violation parameters fieff-, ^CP across the B^ K~^K~K^ Dalitz plot, and 
remove an ambiguity in the solution for the mixing angle [3eff ^ 7r/2 — (ieff^ present in previous 
measurements of sin(2/3ejff) in penguin decays. In our analysis the reflection is removed due to 
interference between CP-even and CP-odd decays that give rise to a cos(2/?e^) term, in addition to 
the sin(2/3ej(f ) terms that come from the interference of decays with and without mixing. Fit results 
for CP parameters are listed in Table 3 for both the K~ K^{'k~^ 'k~) submode and the entire 
sample. At projection plots for signal-weighted events [11], shown in the right plots of Fig. 9, 
clearly show a large phase difference between B^ and B^ decays. The correlation coefficients r 
between the CP parameters are given in Table 3. The global correlation coefficients in the fit to 
the combined sample are 0.02 and 0.08 for j3f.ff and Acp-, respectively. 
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IN dlliC 


i^+i^-K0(7r+7r-) Combined 


Acp 
r 


-0.100 ±0.089 -0.034 ±0.079 ±0.025 
0.354 ±0.083 0.361 ± 0.079 ± 0.037 
0.003 0.013 



Table 3: Average CP-asymmetry parameters for the B K^K^K Dalitz plot. For the combined 
K^K~K^ sample, the first error is statistical and the second is systematic. For the K'^ K~ 
submode, only statistical errors are shown. 




Figure 10: K^K K^{t:^t: ): Change in the Alog(>C) value as a function of fieff- 



Using the K^^-k^ tt~) subsample, we estimate the significance of the nominal result for 

[5 eff compared to the trigonometric reflection where (3 eff — > 7r/2 — f3eff- In a collection of fits with 
both isobar coefficients and CP-asymmetry parameters allowed to vary, we randomize the initial 
parameter values and evaluate the likelihood separation between these two solutions. We find 
A log(£) = 10.4, which excludes the reflection at a significance of 4.6 standard deviations. Note 
that a reflection f3eff — > (3eff ± still remains since we measure the total phase difference between 

and B^ decays {213 eff)- A scan of the change in likelihood as a function of (3 eff is shown in 
Figure 10. 

5 SYSTEMATIC STUDIES 

We study systematic effects on the CP-asymmetry parameters due to fixed parameters in the event- 
selection (rn-Es and IS.E) PDFs. We assign systematic errors by comparing the fit with nominal 
parameters and with parameters smeared by their error, and assign the average difference as the 
systematic error. We account for a potential fit bias using values observed in studies with MC 
samples generated with the nominal Dalitz plot model. We take the largest values of the bias 
observed in these studies as the systematic error. We account for fixed resolution parameters, 
B^ lifetime, B^-B^ mixing and flavor tagging parameters. We also assign an error due to interference 
between the CKM-suppressed b ucd and the favored b — > cud amplitude for some tag-side B 
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decays [26]. Smaller errors due to beam-spot position uncertainty, detector alignment, and the 
boost correction are based on studies done in charmonium decays. In all fits we assume no direct 
CP violation in decays dominated by the 6 — c transition [xcqK^-, D(^g-^K). 



Parameter 


^CP Peff ^CP Peff 


K+K-K^ 

^CP Peff 


Event selection 
Fit Bias 
At, vertexing 
Tagging 
Dalitz model 


0.00 0.01 0.00 0.00 
0.04 0.01 0.04 0.02 
0.02 0.03 0.01 0.01 
0.01 0.00 0.01 0.00 
0.09 0.03 0.09 0.03 


0.003 0.002 
0.004 0.010 
0.010 0.010 
0.021 0.002 
0.011 0.035 


Total 


0.10 0.05 0.10 0.04 


0.025 0.037 



Table 4: Summary of systematic errors on CP-asymmetry parameters. Errors for cpK and foK 
CP-parameters are based on the low-iC+iC~-mass sample. The K~^K~K^ column refers to errors 
on average CP parameters across the Dalitz plot. 

In the fit to low K'^ K~ masses, we extract CP asymmetry parameters for (f)K^ and foK^ 
decays together with isobar parameters for the (^-fC°, while all other isobar parameters are fixed. We 
evaluate the impact of the fixed parameters on the CP violation results using samples of simulated 
events. We compare the reference fit result with the fit that has fixed Dalitz parameters smeared by 
their errors which are taken from Table 1. The average difference in the CP violation parameters 
is taken as the systematic error. We also assign an error due to uncertainty in the resonant and 
non-resonant line-shape parameters. For resonant components this includes the unccrtainity in 
the mass and width of the Xo(1550). For that resonance, we replace our nominal parameters 
with those found by different measurements: = 1.491 GeV/c^, V = 0.145 GeV [17], and 
rrir = 1.507 GeV/c^, F = 0.109 GeV [2]. We take the largest observed difference from the reference 
fit as the systematic error. The non-resonant distributions arc not motivated by theory so we try 
several alternative non-resonant models which omit some of the dependences on K~^Kg and K~Kg 
masses (see Eq. 9): e~""^i2, e~""*i2 + C23e^'f'^^ e~°''^23 ^ and e~"'"i2 + ci^e^'^^^ e~""^i3 . We also study 
the effect of the uncertainty of the shape parameter a on the CP parameters. The non-resonant 
events contribute to the background under the (j) but their shape is determined from the high-mass 
region. We therefore omit the non-rcsonant terms and re-do the low-mass fits, and take the largest 
difference from the reference fit as a systematic error. In the full Dalitz plot fit, the non-resonant 
term is the dominant contribution to the sample and its omission is not reasonable. The mass 
resolution is neglected in the reference fit since it is small compared to the resonant width for all 
Dalitz plot components. We evaluate a potential bias on CP-asymmetry parameters by repeating 
the fit with the Dalitz plot PDF convolved with the mass resolution function and take the difference 
in the CP parameters into the systematic error. 

In the fit to the charge asymmetry in — > (f)K~^ we consider systematic errors due to charge 
asymmetries in tracking and particle identification (0.011), uncertainties in the parameterization 
of the signal Fisher PDF (0.003) and B background content (0.012). We add these contributions 
in quadrature to obtain the total systematic uncertainty on the direct CP violation. 
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6 CONCLUSIONS 



In a fit to K^{7t^tt~) decays, we analyze the Dalitz plot distribution and measure 

the fractions to intermediate states, given in Table 1. Subsequently, we extract CP-asymmetry 
parameters from simultaneous fits to K'^K'K^ final states with the neutral kaon reconstructed 
as Kg — >■ 7r+7r~,K° — 7r°7r°, or K^. We further analyze the K'^K~ phase-space by computing 
moments of Legendre polynomials in (/lET"*" and K~^K~ Kg{7r~^TT~) decays. We find the P-wave 
fraction to be 0.29 it 0.03 averaged over the Dalitz plot, and 0.89 ±0.01 over the ^(1020) resonance 
region (1.0045 < uik+k- < 1-0345 GeV/c^). 

Prom a fit to events at low K+K' masses, we find P^ff = 0.06 ± 0.16 ± 0.05 for B° ^K^ 
and 0.18 it 0.19 it 0.04 for foK^, consistent with our previous measurements [9] and with an 

update of the previous method to the present dataset. We do not observe any significant deviation 
of CP parameters from the Standard Model values P ~ 0.37, Acp = 0. 

In a fit to the full Dalitz plot, we find the CKM angle f3eff = 0.361 ± 0.079 ± 0.037 to be 
compatible with the SM expectation. Additionally, we resolve the trigonometric ambiguity in the 
measurement of Peff at 4.6 standard deviations, which is the first such measurement in penguin 
decays. 
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